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In  the  present  work,  microstructural  and  transport  properties  of  a  three-dimensional  (3D)  microstructure 
of  lanthanum  strontium  manganite  (LSM)  are  deduced  using  dual-beam  focused  ion  beam-scanning  elec¬ 
tron  microscopy  (FIB-SEM)  facility.  A  series  of  two-dimensional  (2D)  cross-sectional  images  are  collected 
from  the  LSM  sample  using  FIB-SEM  and  then  reconstructed  to  3D  structures  from  the  2D  images  in  a 
systematic  approach.  For  the  first  time,  the  effect  of  different  image  processing  steps  including  thresh¬ 
old  value,  median  filter  radius,  morphological  operators,  surface  triangulation,  smoothing  filter,  etc.,  on 
porosity,  internal  surface  area,  electronic  conductivity  and  diffusivity  are  studied.  Variation  of  33%  and 
25%  on  porosity  s  and  internal  surface  area  S,  respectively  is  observed  because  of  improper  selection  of 
threshold  value,  median  filter  radius,  and  morphological  operator.  The  number  of  triangular  surfaces  used 
in  3D  reconstructions  also  varied  the  porosity  s  and  internal  surface  area  S  by  1 4.5%  and  4.4%,  respectively. 

Computational  domains  for  calculating  effective  transport  properties  are  generated  using  body-fitted 
cut-cell  based  finite  volume  meshes  on  reconstructed  3D  volumes.  The  normalized  effective  transport 
properties  are  computed  on  computational  domains  reconstructed  by  the  FIB-SEM  as  well  as  by  a  numer¬ 
ical  model.  For  the  FIB-SEM  reconstruction  case,  the  normalized  effective  properties  in  z-direction  are 
25-44%  smaller  than  those  properties  in  x  and  y  directions.  This  difference  is  significant  and  reveals  the 
anisotropy  in  FIB-SEM  reconstructed  volume  compared  to  numerically  reconstructed  volume.  The  pres¬ 
ence  of  large  crater,  milling  direction  and  smaller  3D  FIB-SEM  reconstructed  volume  could  be  the  main 
reasons  for  this  local  anisotropy. 

©  2010  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

There  has  been  an  increasing  interest  in  understanding  three- 
dimensional  (3D)  microstructures  of  porous  materials  such  as  those 
used  for  electrode  materials  in  solid-oxide  fuel-cells  (SOFC)  [1-11] 
so  as  to  establish  structure-property  and  structure-performance 
relationships.  The  accurate  determination  of  effective  transport 
properties  such  as  gas  diffusivity,  electronic  conductivity,  and 
permeability  is  required  for  designing  highly  efficient  electrodes 
for  energy  conversion  devices  such  as  SOFCs  [2-5].  Microstruc- 
tural  parameters  like  fluid-solid  interfacial  area,  particle  size 
distribution,  pore  size  distribution  as  well  as  local  porosity  or 
material  composition  are  also  important  in  this  regard  [6-9]. 
Hence,  extracting  the  desired  microstructural  information  from  the 
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porous  materials  has  significant  relevance  to  the  fuel  cell  indus¬ 
try. 

Scanning  electron  microscopy  (SEM)  is  a  useful  technique  for 
extracting  two-dimensional  (2D)  images  of  the  microstructures 
but  it  does  not  provide  the  third  spatial  component  (i.e.,  thick¬ 
ness)  of  the  sample  which  is  important  to  find  interconnected 
regions  and  pore  volumes,  shapes  and  sizes.  Typically,  stereolog- 
ical  analysis  is  used  to  obtain  the  3D  microstructures  and  their 
parameters.  In  stereological  analysis,  the  feature  geometry  must  be 
assumed  to  be  that  of  plates,  cylinders,  spheres,  oblate  and  prolate 
spheroids,  or  other  structures  according  to  geometric  probability 
[1  ].  At  the  same  time,  it  is  difficult  to  determine  these  microstruc¬ 
tural  parameters  and  measure  effective  transport  properties  of  a 
porous  electrode  material  experimentally.  Therefore,  an  alterna¬ 
tive  approach  has  been  used  to  construct  such  structures  virtually 
and  to  subsequently  determine  the  microstructural  and  transport 
properties.  Choi  et  al.  [2,1 2,1 3 ],  Kenney  et  al.  [3 ],  and  Metcalfe  et  al. 
[14]  numerically  constructed  the  3D  porous  electrode  structures 
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by  using  measurable  starting  parameters  such  as  the  particle  size 
distribution.  The  advantage  of  such  numerical  models  is  that  they 
can  quickly  generate  a  desired  3D  porous  geometry  from  which 
important  microstructural  parameters  can  be  estimated.  However, 
these  structures  are  based  on  building  the  electrode  geometry  using 
numerical  algorithms  considering  random  placement  of  spheres 
with  some  allowance  of  particle  sintering  through  overlaps.  There¬ 
fore,  these  methods  are  only  approximations  of  the  actual  physical 
microstructure  of  the  porous  electrodes.  Real  microstructures  need 
to  be  produced  from  the  porous  materials  like  SOFC  electrodes  to 
compare  against  the  numerical  ones  and  for  a  better  understanding 
of  the  microstructural  details. 

Tomography  is  one  method  to  generate  the  realistic  3D  images 
from  the  2D  sectional  images  [15].  There  has  been  increasing 
application  of  this  technique  for  microstructure  characterization 
of  fuel  cell  electrodes  using  nano  X-ray  computed  tomography 
(nano-XCT)  [6,16-28]  and  dual  beam  focused  ion  beam-scanning 
electron  microscopy  (FIB-SEM)  [4,5,8,9,29-34].  Nano-XCT  is  a  non¬ 
destructive  technique  used  to  characterize  the  fuel  cell  electrodes, 
semiconductors,  rock  samples,  biological  samples,  etc.  [21].  Izzo 
et  al.  [17-19,23]  used  nano-XCT  to  acquire  series  of  images  at  spa¬ 
tial  resolution  of  42.7  nm  for  3D  reconstruction  of  SOFC  sample 
with  solid  YSZ  electrolyte  and  porous  nickel-yttria  stabilized  zir- 
conia  (Ni-YSZ)  anode  layers.  In  their  work,  porosity  was  calculated 
using  voxel  count  method  and  tortuosity  was  evaluated  by  solving 
Laplace  equation.  They  validated  volume  fraction  calculations  with 
mercury  intrusion  porosimetry.  Lattice  Boltzmann  Method  (LBM) 
was  also  used  to  estimate  the  depletion  of  fuel  gas  in  the  fuel  cell 
electrode.  Grew  et  al.  [24]  have  also  investigated  the  structure  of 
Ni-YSZ  anode  using  nano-XCT  with  a  spatial  resolution  of  38.5  nm. 
They  extracted  tortuosity,  contiguity  and  pore  sizes  and  then  com¬ 
pared  their  results  with  LBM  of  Izzo  et  al.  [18].  Shearing  et  al.  [6] 
have  also  reported  the  nano-XCT  study  to  understand  the  struc¬ 
ture  of  the  Ni-YSZ  electrode  material.  They  performed  nano-XCT  at 
varying  length  scales  of  Ni-YSZ  sample  and  compared  their  struc¬ 
tural  details  with  FIB  tomography.  Gaun  et  al.  [27,28]  studied  the 
structure  of  LSM-YSZ  and  Ni-YSZ  SOFC  composite  electrodes  using 
nano-XCT.  The  results  obtained  in  their  work  were  used  to  quantify 
structural  parameters  like  volume  ratio  of  three  phases,  connected 
porosity,  surface  area  of  each  phase,  etc.  However,  this  nano-XCT 
technique  has  a  significant  short-coming  in  that  a  very  small  sample 
is  required  to  obtain  high  resolution  images,  which  is  often  in  the 
order  of  millimetre  dimensions  or  even  smaller.  Such  limitations  in 
sample  size  can  be  overcome  by  using  a  FIB  for  the  sample  prepara¬ 
tion  [6].  Ostadi  et  al.  [10]  studied  the  3D  reconstructed  images  for 
gas  diffusion  layer  (GDL)  and  a  microporous  layer  (MPL),  relevant 
to  Proton  Exchange  Membrane  (PEM)  fuel  cells.  They  created  a  3D 
image  of  the  GDL-MPL  assembly  through  nano-XCT  with  a  680  nm 
pixel  resolution  for  the  GDL  and  FIB-SEM  with  14nm  pixel  res¬ 
olution  for  the  MPL.  They  determined  the  critical  nano-structural 
features  like  porosity,  characteristic  lengths  and  3D  pore  size  distri¬ 
bution  directly  from  3D  images.  They  also  computed  the  tortuosity 
and  permeability  of  the  MPL  using  LBM. 

FIB,  a  destructive  imaging  technique,  first  came  into  existence 
two  decades  ago  for  nano-machining,  depositing  and  imaging  semi¬ 
conductor  materials.  It  is  a  well-established  technique  to  generate 
3D  microstructures  of  porous  material  with  resolutions  of  nanome¬ 
tres  and  without  any  dependence  on  sample  size  [8,9,35].  FIB 
tomography  is  very  useful  for  studying  SOFC  microstructures,  since 
most  SOFC  electrodes  have  nano-scale  pores.  Wilson  et  al.  [4, 5, 8, 9] 
have  demonstrated  FIB-SEM  by  studying  the  microstructure  of 
composite  electrodes  such  as  Ni-YSZ  cermet  anodes,  LSM-YSZ 
cathodes,  and  (La,Sr)Co03  cathodes  using  dual-beam  FIB-SEM  facil¬ 
ity.  To  aid  in  the  analysis  of  composite  electrodes,  Wilson  et  al.  [4,5] 
used  an  impregnated  sample  with  an  epoxy  to  discriminate  differ¬ 
ent  phases  in  the  images  and  used  inlens  detectors  for  capturing 


back-scattered  electrons  in  the  images.  Wilson  et  al.  [8]  extracted 
the  macroscopic  electrode  parameters  such  as  interface  areas  and 
triple-phase  boundaries  (TPB)  from  the  3D  microstructural  data. 
They  also  discussed  the  calculation  of  critical  microstructural  fea¬ 
tures  such  as  volume  fractions  and  surface  areas  of  specific  phases, 
TPB,  and  the  connectivity  and  tortuosity  of  specific  subphases 
within  SOFC  anode  (Ni-YSZ)  [9].  In  a  similar  study,  Wilson  et  al.  [5] 
reported  the  3D  characterization  and  quantification  of  important 
structural  parameters  from  LSM-YSZ  cathode  and  identification 
of  LSM-YSZ-pore  three-phase  boundaries.  They  also  studied  the 
effect  of  composition  of  cathodes  on  performance  of  SOFC  by  cor¬ 
relating  the  3D  microstructure  and  polarization  resistance  [4].  In 
the  present  study,  similar  to  Wilson’s  set-up,  the  FIB-SEM  is  used  to 
record  the  2D  sectional  images  of  lanthanum  strontium  manganite 
(LSM),  a  typical  cathode  material  of  SOFC. 

Following  the  FIB-SEM  work  reported  by  Wilson  et  al.,  there  has 
been  a  growing  interest  on  characterization  of  the  microstructure  of 
electrode  materials  by  various  techniques  and  subsequent  extrac¬ 
tion  of  micro-structural  property  information.  Kenney  et  al.  [3] 
computed  the  TPB  length,  surface  area  and  pore  size  from  numer¬ 
ical  reconstruction  of  composite  SOFC  electrodes.  They  found  that 
the  TPB  length  evaluated  for  mono-sized  particles  with  same  mean 
particle  sizes  was  as  much  as  40%  higher  compared  to  real  parti¬ 
cle  size  distributions.  Iwai  et  al.  [1 1  ]  also  evaluated  the  TPB  length 
density  and  tortuosity  factors  from  SOFC  anode  (Ni-YSZ)  using 
the  dual-beam  FIB-SEM  equipped  with  an  energy  dispersive  X- 
ray  spectroscopy  (EDX)  unit.  They  estimated  the  TPB  density  using 
the  centroid  and  volume  expansion  methods  and  tortuosity  factors 
with  lattice  Boltzmann  and  random  walk  methods.  They  observed 
the  excellent  agreement  between  the  two  methods  used  for  calcu¬ 
lating  each  parameter. 

As  elucidated  in  this  previous  paragraph,  there  is  a  significant 
interest  by  different  research  groups  to  compute  micro-structural 
properties  from  the  images  derived  from  FIB-SEM  studies.  Notwith¬ 
standing  the  challenges  in  fine-tuning  of  the  experimental 
parameters  to  generate  high  quality  2D  data,  the  procedure  for 
the  post-processing  of  images  to  generate  representative  3D  struc¬ 
ture  and  the  methodology  for  subsequent  generation  of  a  mesh  for 
the  computation  of  transport  properties  affect  the  final  outcome. 
Processing  of  images  is  required  before  stacking  the  2D  images  to 
obtain  a  3D  reconstructed  image  and  obtaining  the  microstructural 
parameters  like  porosity,  composition,  feature  size,  shape  and  con¬ 
nectivity.  The  processing  steps  like  image  registration,  algebraic 
operations,  morphological  operations,  filtering  and  segmentation, 
etc.,  are  required  prior  to  3D  reconstruction  to  align  the  images;  to 
remove  any  noise  and  to  label  the  different  phases.  Recently,  Joos 
et  al.  [34]  described  the  effect  of  image  processing  steps  on  the  cal¬ 
culation  of  surface  area,  volume/porosity  fraction  and  tortuosity. 
Their  work  is  mainly  focused  on  the  selection  of  proper  threshold 
value  to  differentiate  the  pore  space  and  solid  matrix.  In  the  present 
work,  we  examine  the  effect  of  different  image  processing  steps  in 
calculating  porosity,  internal  surface  area,  diffusivity,  etc.  A  number 
of  statistical  [36-40]  and  process-based  [38,41]  models  are  devel¬ 
oped  for  reconstructing  3D  images  from  2D  thin  section  images. 
In  the  present  study,  a  modified  and  improved  marching  cubes 
algorithm  [42-46]  is  used  for  generating  3D  reconstructed  images. 
However  for  calculating  effective  transport  properties,  we  con¬ 
verted  3D  reconstructed  microstructures  into  body-fitted  cut-cell 
based  finite-volume  meshes  using  MicroFOAM  [12],  which  is  a  suit¬ 
ably  modified  version  of  the  open-source  CFD  code  OpenFOAM® 
[47]. 

The  existing  body  of  literature  on  extraction  of  microstruc¬ 
ture  of  porous  material,  like  LSM,  using  FIB-SEM  seldom  describes 
the  impact  of  various  image  processing  steps  and  algorithm  on 
the  phase  differentiation,  in  particular  the  porosity  and  internal 
surface  areas.  The  main  focus  of  this  paper  is  to  describe  a  system- 
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atic  approach  in  quantifying  the  processed  images  from  FIB-SEM 
and  subsequently  computing  the  transport  properties  of  porous 
electrode,  like  electronic  conductivity  and  diffusivity.  Further,  a 
systematic  analysis  is  reported  for  the  first  time  in  this  work,  where 
the  effects  of  various  image  processing  steps  including  threshold 
value,  median  filter  radius,  morphological  operators,  surface  trian¬ 
gulation,  smoothing  filter,  etc.,  are  demonstrated  on  the  computed 
porosity,  internal  surface  area,  electronic  conductivity  and  diffu¬ 
sivity.  A  sensitivity  study  is  also  performed  where  the  computed 
transport  properties  are  compared  between  those  obtained  from 
reconstructed  FIB-SEM  images  and  the  numerical  reconstructed 
images  assuming  random  structures. 

2.  Materials  and  methods 

2.1.  Sample  preparation 

The  porous  cathode  was  deposited  using  (Lao.8Sr0.2)o.95Mn03±(5 
powder  (LSM,  Praxair)  onto  a  thick  (Y203  )o.o8(Zr02  )o.92  (YSZ,  Tosoh) 
electrolyte  disk.  The  LSM  powder  particle  size  was  determined  from 
a  light  scattering  technique  to  have  a  mean  diameter  of  0.73  |xm 
with  a  standard  deviation  of  0.3  [xm.  First,  the  YSZ  powder  was 
pressed  in  a  steel  die  set  using  a  Carver  10  ton  press  at  a  pressure 
of  1 5,000  psi.  The  electrolyte  disk  was  then  sintered  to  1400  °C  for 
5h.  Once  cooled,  the  electrolyte  surface  was  polished  with  240, 
400  and  600  grit  silicon  carbide  paper  and  re-fired  to  1400  °C.  After 
the  sintering  processes,  the  dimensions  of  the  electrolyte  were 
found  to  be  approximately  12.9  mm  in  diameter  and  0.7  mm  thick. 
Next,  LSM  powder  was  mixed  with  ethanol,  poly(vinylpyrrolidone) 
and  5  mm  diameter  YSZ  milling  media  and  the  mixture  was  ball 
milled  for  20  h.  The  resulting  LSM  ink  was  then  pulse  sprayed 
using  an  automated  air-brush  onto  the  YSZ  disk.  The  combined 
electrode/electrolyte  was  then  sintered  to  1 150  °C  for  2  h. 

The  porosity  of  electrodes  fabricated  using  this  method  was 
estimated  from: 


where  peff  is  the  effective  density  of  the  electrode,  estimated  based 
on  the  volume  of  the  electrode  and  the  deposit  weight,  and  pbulk 
is  the  bulk  density  of  an  LSM  particle,  equivalent  to  the  electrode 
density  for  an  electrode  with  zero  porosity.  The  thickness  of  the  sin¬ 
tered  electrode  was  16  |xm,  measured  from  an  SEM  cross-section 
image,  the  deposit  weight  of  the  electrode  was  7.1  mg  and  the 
bulk  density  was  6.39 gem-3.  With  these  properties,  the  average 
porosity  of  the  electrode  was  estimated  to  be  42.5%. 

2.2.  Acquisition  of  2D  images  from  FIB-SEM 

Flere  we  describe  the  steps  followed  to  obtain  FIB-SEM  images. 
First,  a  layer  of  10  nm  chromium  and  40  nm  gold  was  sequentially 
deposited  onto  the  LSM  electrode  sample  using  a  planar  mag¬ 
netron  sputtering  system.  Then,  the  top  surface  of  the  sample  was 
inspected  using  a  Zeiss  N-Vision  40  crossbeam  workstation.  The 
N-Vision  40  has  the  capability  of  doing  both  FIB  milling  and  SEM 
imaging,  and  is  therefore  known  as  a  dual-beam  FIB-SEM  facility. 
Next,  a  trapezoidal  trench  of  width  19.97  [xm,  height  34.09  [xm, 
depth  25  |xm  and  base  angle  30°  was  milled  into  the  LSM  sam¬ 
ple  under  an  accelerated  voltage  of  30  kV  and  a  milling  current  of 
27  nA.  This  trench  helps  in  better  exposure  of  the  section  of  inter¬ 
est  and  contrast  adjustment  purpose  at  the  time  of  taking  series 
of  cross-sectional  SEM  images.  Fig.  1(a)  shows  the  location  of  the 
section  considered  on  the  LSM  sample  and  Fig.  1(b)  shows  low- 
magnification  image  of  LSM  sample  with  a  trench.  Thereafter,  a 
layer  of  1  |xm  platinum  was  deposited  on  the  LSM  sample  sur¬ 
face  for  a  rectangular  area  of  11.94|xm  width  and  12  |xm  height 


as  shown  in  Fig.  1(b)  under  an  accelerated  voltage  of  30  kV  and  a 
milling  current  of  300  pA  using  deposition  mode  of  the  N-Vision 
40  system.  The  purpose  of  this  platinum  layer  was  to  correct  the 
tilt  and  alignment  while  taking  the  serial  sectioning  images  into 
consideration.  Next,  the  electron  beam  was  used  to  image  the  side 
wall  of  the  trench  (platinum  coated  side)  at  an  angle  of  54°  from 
the  parallel  to  the  side  wall.  A  trapezoidal  area  of  width  8.99  [xm, 
height  8.99  |xm  and  base  angle  20°  was  selected  on  the  platinum 
coated  surface  for  milling  the  surface  for  taking  the  serial  section¬ 
ing  images.  Then  a  40  nm  thick  layer  was  milled  in  the  selected 
area  under  an  accelerated  voltage  of  30  kV  and  a  milling  current 
of  700  pA.  A  series  of  2D  images  of  10.94nm  resolution  were  col¬ 
lected  from  the  LSM  sample  by  milling  the  sample  layer  by  layer. 
A  set  of  157  images  were  collected.  Fig.  1(c)  shows  the  top  view 
of  the  trench  and  milled  volume  from  low  resolution  to  high  res¬ 
olution.  The  bottom  surface  of  the  milled  volume  was  found  to 
be  not  perfectly  smooth,  as  observed  from  the  enlarged  view  in 
Fig.  1(c).  Fig.  1(d)  depicts  the  trench  and  milled  volume  at  an  angle 
of  55°  from  the  sample  table  to  FIB  milling  direction.  The  serial  sec¬ 
tion  of  slice  number  151  is  shown  in  the  enlarged  figure,  in  which 
grey  colour  indicates  the  LSM  solid  matrix  and  black  colour  indi¬ 
cates  the  pore  space.  The  images  shown  in  Fig.  1  were  obtained 
using  the  secondary  electron  signal  detector.  Imaging  with  the  elec¬ 
tron  beam  was  done  at  voltages  around  2-5  kV  whereas  milling  of 
trenches  was  performed  with  Ga  ion  beam  operating  up  to  30  kV 
and  a  beam  current  of  27  nA.  These  beam  currents  and  voltages  are 
optimized  for  LSM  sample  for  minimizing  the  charging  effect  and 
electron/ion  beam  damage  to  surface.  The  penetration  of  Ga  ions 
can  be  checked  using  open  source  software:  SRIM  (The  Stopping 
and  Range  of  Ions  in  Matter).  If  the  ion  beam  operates  at  voltage  and 
current  greater  than  30  kV  and  27  nA,  the  sample  becomes  amor¬ 
phous  due  to  the  penetration  of  Ga  ions.  Ga  ions  may  penetrate 
up  to  20-40  nm  within  the  substrate.  This  is  more  of  a  concern  for 
really  thin  samples  (TEM  work)  [1  ]. 

2.3.  3D  stacking  or  reconstruction 

A  3D  reconstruction  of  the  porous  LSM  structure  was  obtained 
by  stacking  the  2D  cross  sectional  SEM  images  in  3D  space.  In  this 
work,  most  of  the  image  processing  and  3D  reconstruction  was 
carried  out  using  the  Avizo  Fire  edition  version  6.2  and  XSkele- 
ton  extensions  (VSG,  Visualization  Sciences  Group,  Inc.,  Burlington, 
MA).  Pre-processing  of  images  was  required  to  get  a  good  qual¬ 
ity  3D  reconstruction.  The  first  step  in  creating  3D  microstructure 
involves  alignment  of  collected  serial  sectioned  2D  images  using 
some  fiduciary  markers  as  objects  of  reference.  Fiduciary  marking 
(reference  points)  used  for  alignment  of  images  is  shown  in  Fig.  1  (e). 
Four  points  were  located  outside  the  milled  volume  on  the  images 
as  fiduciary  markers  and  these  points  were  then  adjusted  to  place 
all  images  at  the  common  point  of  observation.  After  executing 
this  fiduciary  marking  alignment,  least-square  algorithm  is  used 
once  again  for  alignment  to  remove  any  errors  in  the  alignment 
based  on  grey  values.  Proper  alignment  of  images  without  down 
sampling  was  required  while  stacking  the  images  to  get  accurate 
3D  reconstruction.  A  volume  of  9.12  |xm  x  4.75  |xm  x  4.08  [xm  was 
selected  from  the  3D  data.  The  selected  volume  should  be  avail¬ 
able  in  every  aligned  images  starting  from  the  first  image  to  last 
sectioned  image.  Since  the  resolution  and  contrast  in  the  sectioned 
images  were  reduced  progressively  from  the  first  to  the  last  section, 
a  process  for  changing  the  range  of  pixel  intensity  values,  i.e.,  nor¬ 
malization  process  was  carried  out.  The  normalized  image  of  slice 
51  is  shown  in  Fig.  2(a).  Then,  the  normalized  image  was  processed 
to  find  different  phases  of  the  sample  by  segmentation  process. 
The  segmentation  process  consists  of  separating  and  identifying 
the  different  segments  of  a  digital  image  based  on  the  intensity  of 
its  pixels.  This  step  is  carried  out  by  adjusting  the  threshold  of  the 
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Fig.  1.  Low-magnification  image  of  LSM  sample:  (a)  location  of  section  considered  in  the  LSM  sample;  (b)  trench  in  the  sample  is  shown  at  stage  angle  of  54°  along  with 
the  region  showing  the  platinum  coated  surface,  this  portion  of  this  coated  surface  is  taken  for  milling  the  material  for  3D  reconstruction;  (c)  the  top  view  of  the  trench 
along  with  the  volume  taken  for  3D  reconstruction;  enlarged  view  in  the  left  shows  the  high  resolution  image  of  the  same  area  and  the  enlarged  view  in  right  shows  the 
bottom  surface  of  region  taken  for  reconstruction  after  milling;  (d)  the  cross-sectional  view  of  the  milled  surface  taken  for  reconstruction;  and  (e)  the  fiduciary  marks  used 
for  alignment  of  series  of  images. 


pixel  intensity  values.  In  this  set  of  images,  the  threshold  value  was 
kept  as  60,  i.e.,  pixel  intensity  value.  The  pixels  that  are  below  60  in 
the  image  were  kept  as  black  and  the  pixels  which  are  in  the  range 
of  60-255  were  kept  as  white.  The  black  colour  in  the  segmented 
image  as  shown  in  Fig.  2(b)  indicates  the  pore  and  the  white  colour 
indicates  the  solid  matrix  (grain),  that  means  two  phases  in  the 
images  were  labeled  as  pore  (void)  space  and  solid  matrix,  respec¬ 
tively.  Nonlinear  digital  filtering  methods  such  as  3D  median  filter 
was  then  used  to  reduce  the  noise  in  images.  This  process  improved 


the  edge  detection  on  all  images  which  is  helpful  for  the  subsequent 
processing  stages.  Further,  the  smoothing  of  images  was  carried 
out  to  reduce  the  noise  within  an  image  and  thus  produce  a  less 
pixelated  image.  Fig.  2  shows  the  change  in  the  image  of  slice  num¬ 
ber  51  after  processing,  as  described  here.  Comparing  the  bottom 
left  corner  of  normalized  image  (Fig.  2(a))  and  the  image  result¬ 
ing  after  smoothening  steps  described  here  (Fig.  2(e)),  it  appears 
that  certain  portion  of  the  solid  matrix  is  incorrectly  assigned  as  a 
pore  space.  Due  to  presence  of  dual  beam  imaging  artifacts,  such  as 
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Fig.  2.  Steps  used  in  preprocessing  the  images  before  reconstruction  is  shown  on  cross  sectional  image  of  slice  number  51 :  (a)  grey  scale  image  after  normalization  process 
which  has  been  corrected  for  shadowing,  where  pore  space  is  dark  and  matrix  is  bright;  (b)  binary  image  computed  by  thresholding  the  image  in  A,  where  pore  space  is  black 
and  matrix  is  white;  (c)  binary  image  computed  by  applying  median  filter  to  image  in  B,  where  pore  space  is  white  and  matrix  is  blue;  (d)  surface  cut  of  3D  reconstructed 
image  after  stacking  using  marching  cubes  algorithm,  the  blue  colour  indicates  the  solid  matrix  and  transparent  colour  indicates  the  pore  or  void  space;  (e)  surface  cut  of  3D 
reconstructed  image  after  applying  smoothing  filter  to  remove  the  noise  in  the  image  and  edges;  and  (f)  the  difference  in  the  pore  and  solid  matrix  phases  before  applying 
smoothing  filter  and  after  smoothing  filter  on  3D  reconstructed  image.  (For  interpretation  of  the  references  to  colour  in  this  figure  legend,  the  reader  is  referred  to  the  web 
version  of  the  article.) 


curtaining,  detector  noise,  the  normalized  image  does  not  clearly 
represent  the  actual  surface  morphology  of  the  slice.  It  is  to  be 
noted  that  in  the  present  study,  the  sample  was  not  infiltrated  by 
epoxy.  In  some  cases,  infiltration  is  a  better  option  because  it  can 
give  better  contrast  and  prevents  the  solid  phase,  which  may  be 
out  of  the  measurement  plane,  to  appear  in  the  images.  Infiltration 
does  not  change  the  segmentation  needs,  it  merely  brings  in  more 
“greylevel  resolution”  between  pore  and  solid  (electrode)  phase. 
The  method  of  visual  inspection  used  in  our  study  yields  a  satisfac¬ 
tory  segmentation  so  that  infiltration  was  not  necessary  in  our  case. 
Segmentation  methods  such  as  level  set  method  [48]  might  provide 
more  accurate  segmentation  but  at  the  expense  of  a  considerable 
increase  in  computational  requirements.  Avizo  also  provides  such 
method  with  watershed  segmentation  (which  is  usually  applied  to 
the  image  gradient)  as  a  robust  and  community-accepted  alterna¬ 
tive.  Since  the  sample  has  only  two  phases,  segmentation  based  on 
threshold  value  is  more  appropriate  and  easy  method  to  apply.  In 
addition,  we  applied  median  filter  and  smoothing  filter  for  satisfac¬ 
tory  segmentation  to  remove  segmentation  errors. 

Stacking  of  the  processed  images  can  give  the  3D  reconstructed 
image  of  LSM  sample.  An  improved  marching  cubes  algorithm  was 
used  for  extracting  a  polygonal  mesh  of  an  isosurface  from  a  3D 
voxels  data  [42-46].  This  was  an  inbuilt  algorithm  available  in  the 
Avizo  Fire  edition,  in  which  the  3D  data  was  processed  in  scan¬ 
line  order  and  calculated  triangle  vertices  using  linear  interpolation 
to  create  triangle  models  of  constant  density  surfaces,  i.e.,  polygo¬ 
nal  isosurfaces  [42].  Fig.  3  represents  the  cross  sectional  view  of 
3D  reconstructed  sample  in  three  directions.  The  coloured  part 
shows  the  LSM  solid  matrix  and  transparent  part  represents  the 
pore  space. 


Fig.  3.  Three  different  cross  sectional  views  of  the  3D  reconstructed  sample. 


2.4.  Calculation  of  effective  transport  properties 

The  net  transport  of  a  species  within  a  composite  structure, 
although  governed  by  the  intrinsic  transport  coefficient,  is  mod¬ 
ified  due  to  the  tortuous  pathways  and  the  presence  of  blocking 
phases  within  a  porous  volume.  Thus,  the  effective  transport  coef¬ 
ficient  (aeff)  for  a  composite  material  can  be  related  to  the  intrinsic 
transport  coefficient  via  some  function  of  the  volume  fraction  of 
the  pertinent  phase  responsible  for  the  transport.  The  functionality 
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Fig.  4.  Reconstructed  3D  microstructure  of  LSM  sample:  (a)  3D  reconstruction  showing  the  solid  matrix  of  LSM  sample  in  coloured  part  and  other  pore  space  in  transparent: 
(b)  3D  reconstruction  of  pore-network  structure  showing  the  solid  matrix  of  LSM  sample  in  transparent  and  other  pore  space  with  colour;  and  (c)  3D  reconstruction  showing 
the  solid  matrix  of  LSM  sample  in  green  colour  and  pore  space  in  blue  colour.  (For  interpretation  of  the  references  to  colour  in  this  figure  legend,  the  reader  is  referred  to  the 
web  version  of  the  article.) 


depends  on  the  microstructure.  To  accurately  compute  the  effective 
transport  properties  for  our  sample,  the  3D  reconstructed  volumes 
of  the  solid  LSM  phase  and  pore  spaces  were  converted  into  body- 
fitted  cut-cell  based  finite-volume  meshes  using  MicroFOAM  [12], 
which  is  a  suitably  modified  version  of  the  open-source  CFD  code 
OpenFOAM®  [47]. 

Transport  of  the  electrons  and  gas  species  within  the  phase 
domains  of  the  electrode  sample  is  assumed  to  be  dominated  by 
diffusion  process.  For  each  phase,  the  relevant  potential  field  0  is 
governed  by  the  diffusive  processes  as  follows: 

V-(aV0)  =  O,  </>e.f2phase,  (2) 

where  a  denotes  the  bulk  diffusion  coefficient  and  i2phase  is  the 
computational  volume  for  each  phase.  Note  that  each  phase  vol¬ 
ume  (*f2phase)  is  embedded  in  the  electrode  sample  volume  (£2) 
so  that  &  =  £2i SM IJ  ^pore-  Dirichlet  (fixed  value)  boundary  condi¬ 
tions  are  imposed  on  the  top  and  bottom  boundaries  of  each  phase, 
i.e.,  0top  =  1  and  0bot  =  0.  Symmetry  (zero  gradient)  boundary  con¬ 
ditions,  i.e.,  (  30/3n)  =  O,  are  presumed  on  the  sides  of  each  phase 
boundaries  as  well  as  solid  and  pore  phase  interfaces.  To  solve  the 
diffusion  equation  (i.e.,  Laplace  equation)  of  the  applied  potential 
field  (0)  for  each  phase  domain  ,f2pbase,  the  finite  volume  method 
(FVM)  was  used.  From  the  solution  of  the  diffusion  equation  in  Eq. 
(2),  the  effective  transport  coefficient  for  gas  diffusion  (Deff)  and 
the  effective  electronic  conductivity  (/ceff)  can  be  obtained. 

The  local  transport  flux  q  for  each  phase  can  be  related  to  the 
generalized  local  concentration  gradient  such  that 


(3) 


where  the  potential  field  0  is  calculated  from  Eq.  (2).  The  average 
(surface  integrated)  local  flux  across  the  surface  boundary  of  each 
phase  should  be  satisfied  such  that 
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Lienee  the  normalized  effective  transport  coefficient  can  be  eval¬ 
uated  by 


/9<W#/9n)ds  m 

a  (A  0/L)S  ’  1  j 

where  A 0  represents  the  potential  difference,  L  is  the  selected 
height  of  the  computational  domain  and  5  denotes  the  selected 
boundary  surface  area  of  each  phase.  Note  that  ( aeff/a )  =  (Deff/D) 
represents  molecular  diffusion  ratio  for  the  pore  phase  and 
(aeff/a 0  =  (keff//<)  represents  electrical  conductivity  ratio  for  the  LSM 
phase. 


3.  Results  and  discussions 


Fig.  4  shows  the  resulting  3D  reconstructed  image  pro¬ 
duced  after  stacking  the  2D  images  taken  from  dual-beam 
FIB-SEM  system.  Only  a  reduced  box  shape  domain  of 
4.98  |jim  x  3.04  |xm  x  4.08  [xm  is  displayed  for  illustration.  The 
coloured  part  in  Fig.  4(a)  is  the  LSM  solid  matrix  and  the  transpar¬ 
ent  part  is  pore  space.  The  image  illustrates  the  3D  microstructure 
of  LSM  and  pore  spaces,  and  their  local  surface  curvatures.  A  better 
representation  of  the  pore  space  is  shown  in  Fig.  4(b)  where  the 


q  =  aV0, 
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Fig.  5.  Reconstructed  3D  microstructure  of  LSM  sample  with  2D  view  of  its  three  sides. 


surfaces  are  bound  to  the  void  space.  The  coloured  portion  in 
this  figure  represents  the  pore  space  and  the  transparent  part 
represents  the  LSM  solid  matrix.  This  is  the  3D  image  of  an  actual 
pore  space.  Fig.  4(c)  shows  the  3D  reconstructed  image  of  LSM 
cathode  with  different  colours  assigned  to  the  solid  matrix  and 
pores.  The  pore  space  and  solid  matrix  of  the  LSM  sample  are  easily 
distinguishable  in  this  figure.  Fig.  5  depicts  the  3D  reconstruction  of 
the  LSM  sample  with  corresponding  side  views.  Flere  the  coloured 
part  indicates  the  LSM  solid  matrix  and  transparent  part  indicates 
the  pore  space. 

3.1.  Impact  of  post-processing  on  the  determination  of 
microstructural  properties 

The  different  post-processing  steps  described  previously  affect 
the  measurements  of  the  properties  of  the  electrode;  in  the  follow¬ 
ing,  we  focus  on  the  porosity  s  and  the  internal  surface  area  S.  The 
internal  surface  area  is  calculated  by  counting  the  number  of  ele¬ 
ment  faces  at  the  interface  between  the  solid  and  the  pore  phases. 
In  the  case  of  a  Cartesian  mesh,  for  e.g.,  when  a  3D  image  is  gen¬ 
erated  by  stacking  2D  SEM  images  together,  counting  voxel  faces 
introduces  a  systematic  bias.  For  a  single  sphere  for  instance,  the 
surface  area  is  overestimated  by  a  factor  1.5  [13]. 

In  the  present  analysis,  firstly,  the  effect  of  the  threshold  value 
t,  from  a  range  of  pixel  intensity  values  of  images  of  0-255,  used 


to  discriminate  the  solid  phase  from  the  pore  phase  is  presented 
in  Fig.  6(a)  for  the  porosity  and  in  Fig.  6(b)  for  the  internal  sur¬ 
face  area.  Both  measured  quantities  depend  significantly  on  the 
choice  of  the  threshold  value.  Between  t=  57  and  t=58,  the  esti¬ 
mated  porosity  increases  by  approximately  33%.  This  is  due  to  an 
incorrect  choice  of  threshold,  for  which  large  portions  of  the  image 
are  assigned  to  the  wrong  phase.  Visual  inspection  of  a  few  ran¬ 
domly  chosen  slices  is  required  to  determine  the  proper  threshold 
value.  For  the  present  geometry,  visual  inspection  yields  an  optimal 
threshold  value  around  t=59  or  t=60.  Between  these  two  thresh¬ 
old  values,  the  variation  in  porosity  and  internal  surface  area  is  only 
4%. 

A  significant  amount  of  noise  is  present  in  the  binary  images 
that  appears  as  artificial  pore  speckles  inside  the  solid  phase.  Direct 
application  of  a  median  filter  removes  the  noise  but  tends  to  trans¬ 
form  it  into  a  large  artificial  pore  space  (indicated  by  red  in  Fig.  7). 
Applying  a  morphological  opening  (MO)  operation  [49]  to  the 
binary  image  beforehand  significantly  reduces  the  amount  of  such 
noise.  An  morphological  opening  is  used  to  remove  some  of  the 
pixels  from  the  edges  of  regions  of  pore  space.  Fig.  8(a)  and  (b) 
shows  the  effect  of  the  radius  r  of  the  median  filter,  with  and  with¬ 
out  MO,  for  the  porosity  and  the  internal  surface  area,  respectively. 
As  expected,  not  applying  the  MO  operation  significantly  overes¬ 
timates  the  porosity  (by  about  10%  for  r= 2  pixel)  and  the  internal 
surface  area  (by  about  25%  for  r  =  2  pixel).  The  radius  of  the  median 
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Fig.  6.  Variation  of  measured  properties  as  a  function  of  the  threshold  value:  (a)  porosity  and  (b)  internal  surface  area.  The  reference  for  the  calculation  of  the  variations  in 
porosity  and  internal  surface  area  is  based  on  a  threshold  value  of  t=  60,  a  radius  r  =  1  pixel  for  the  median  filter  and  application  of  the  morphological  opening  operation. 


Fig.  7.  Cross  sectional  binary  images  of  slice  number  51 :  (a)  after  applying  a  median 
filter  of  radius  r= 2  pixel;  (b)  after  applying  a  morphological  opening  and  then  a 
median  filter  of  radius  r=  2  pixel.  The  regions  highlighted  in  red  correspond  to  arti¬ 
ficial  pore  regions  generated  by  the  filter  due  to  the  noise  in  the  original  binary 
images.  (For  interpretation  of  the  references  to  colour  in  this  figure  legend,  the  reader 
is  referred  to  the  web  version  of  the  article.) 

filter  does  not  affect  the  measured  porosity  significantly.  However, 
the  internal  surface  area  varies  by  as  much  as  25%  between  a  filter 
with  a  radius  of  1  pixel  and  one  with  a  radius  of  5  pixel. 

During  the  generation  of  the  surface  mesh,  the  measured  geom¬ 
etry  is  interpolated.  The  structure-dependent  properties  of  the 
material  will  be  more  or  less  altered,  depending  on  the  quality  of 
the  interpolation.  The  surface  mesh  (i.e.,  triangulation)  is  generated 
for  each  phase  of  the  3D  image  of  9.12  [xm  x  4.75  [xm  x  4.08  |xm. 
Fig.  9  demonstrates  the  porosity  (%),  internal  surface  area  (|xm2), 
and  internal  surface  area  density  (|xm2  |xm-3)  as  a  function  of  the 
number  of  surface  triangulation  before  and  after  applying  smooth¬ 
ing  filter.  Fig.  9(a)-(c)  clearly  shows  that  as  the  number  of  surface 
triangulation  of  the  surface  mesh  increases,  the  porosity  and  the 
internal  surface  area  tend  towards  asymptotic  values  when  the 
number  of  surface  triangulation  exceeds  500,000.  The  smoothing 
filter  of  the  surface  triangulation  results  in  the  internal  surface  area 
changing  from  14.5%  to  1.1%,  as  shown  in  Fig.  9(b)  and  (c),  and 
the  porosity  from  4.4%  to  0.03%,  as  shown  in  Fig.  9(a).  Since  the 
porosity  affects  the  effective  gas  diffusivities,  any  changes  in  surface 
triangulation  that  affects  porosity  also  affects  the  numerically  com¬ 
puted  effective  diffusivity.  Fig.  10  shows  the  normalized  effective 
gas  diffusivity  as  a  function  of  the  number  of  surface  triangulation 
before  and  after  applying  smoothing  filter.  Fig.  10(a)  -  (c)  shows 
that  the  normalized  effective  gas  diffusivities  in  x,y,  and  z  directions 


Fig.  8.  Variation  of  measured  properties  as  a  function  of  radius  of  the  median  filter, 
with  and  without  morphological  opening:  (a)  porosity  and  (b)  internal  surface  area. 
The  reference  for  the  calculation  of  the  variations  in  porosity  and  internal  surface 
area  is  based  on  a  threshold  value  of  t=  60,  a  radius  r  =  1  pixel  for  the  median  filter 
and  application  of  the  morphological  opening  operation. 


(DJff/D,  Dff/D  and  Dff/D)  are  asymptotically  converging  with  the 
number  of  surface  triangulation  greater  than  500,000. 

Finally,  Fig.  1 1  displays  the  normalized  effective  properties  as 
a  function  of  the  number  of  surface  triangulation  after  applying 
smoothing  filter.  The  normalized  effective  electronic  conductivi¬ 
ties  in  x,  y,  and  z  directions,  i.e.,  kff/k,  kf{ /k  and  are  given 
in  Fig.  11(a)  whereas  the  normalized  effective  gas  diffusivities  in 
x,  y,  and  z  directions,  i.e.,  Dff/D,  Dff/D  and  Dfff/D,  are  shown 
in  Fig.  11(b).  Similar  to  the  geometric  properties,  the  normalized 
effective  transport  properties  are  asymptotically  converging  as  the 
number  of  surface  triangulation  increases.  It  is  noteworthy  that  the 
effective  properties  in  the  z-direction  is  significantly  different  than 
those  inx-  andy-direction  indicating  anisotropy  of  the  overall  sam¬ 
ple  or  at  least  locally  in  the  sample.  This  issue  is  further  discussed 
in  the  Section  3.3. 


3600 


N.S.K.  Gunda  et  al.  /  Journal  of  Power  Sources  196(2011)  3592-3603 


Number  of  surface  triangulation 


Fig.  9.  Geometric  properties  as  a  function  of  the  number  of  surface  triangulation  before  and  after  applying  smoothing  filter:  (a)  porosity  (%);  (b)  internal  surface  area  (pm2); 
and  (c)  internal  surface  area  density  (pm2  pm-3). 
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Fig.  10.  Normalized  effective  gas  diffusivity  as  a  function  of  the  number  of  surface  triangulation  before  and  after  applying  smoothing  filter:  (a)  normalized  effective  gas 
diffusivity  in  x-direction,  Dff/D;  (b)  normalized  effective  gas  diffusivity  iny-direction,  DJff/D;  and  (c)  normalized  effective  gas  diffusivity  in  z-direction,  DfjF/D. 


3.2.  Porosity  and  internal  surface  area 

The  total  volume  (pore  space +  LSM  solid  matrix)  of  the  3D 
reconstructed  structure  is  separated  into  pore  and  LSM  phase 
volumes  through  the  image  post-processing  mentioned  in  the  pre¬ 


vious  section.  Using  the  threshold  value  of  60,  the  porosity  (i.e., 
volume  fraction  of  pore  space)  and  the  volume  fraction  of  LSM 
phase  are  calculated  as  50.6%  and  49.4%,  respectively.  The  porosity 
computed  from  the  FIB-SEM  sample  is  larger  than  the  one  esti¬ 
mated  during  the  fabrication.  This  is  likely  due  to  the  presence  of 
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Fig.  11.  Normalized  effective  properties  in  x,  y,  and  z  directions  as  a  function  of 
the  number  of  surface  triangulation  after  applying  smoothing  filter:  (a)  normalized 
effective  conductivity  in  x,  y,  and  z  directions,  kffz/k  and  (b)  normalized  effective 
gas  diffusivity  in  x,  y,  and  z  directions,  D*ny  Z/D. 


Fig.  12.  Craters  in  the  LSM  sample. 


crater  that  makes  the  sample  more  porous  than  expected  locally. 
Fig.  1 2  shows  the  larger  crater  in  the  2D  cross-sectional  image  (slice 
number  94)  of  LSM.  The  gas-solid  internal  surface  area  per  unit 
volume  (i.e.,  the  internal  surface  area  density,  S/i?)  is  calculated  as 
4.13  |jim2  |jinrr3.  Since  LSM  sample  has  only  two  phases,  the  inter¬ 
nal  surface  area  is  the  interface  between  pore  and  LSM  phases. 


Fig.  13.  3D  computational  domain  considered  for  numerical  simulations 
to  calculate  the  transport  properties:  (a)  ‘Sample  A’  has  the  dimension 
of  9.12  p,m  x  4.75  (Jim  x  4.08  p,m  and  (b)  ‘Sample  B’  has  the  dimension  of 
4.98  |jim  x  3.04  |xm  x  4.08  |jim.  The  green  colour  surface  indicates  LSM  phase  domain 
where  as  the  blue  colour  surface  with  face  grids  denotes  the  pore  phase  domain.  (For 
interpretation  of  the  references  to  colour  in  this  figure  legend,  the  reader  is  referred 
to  the  web  version  of  the  article.) 

The  total  surface  area  densities  (i.e.,  sum  of  internal  and  boundary 
surface  areas  per  unit  volume)  for  pore  and  LSM  phases  are  com¬ 
puted  as  4.91  |jim2  p.m-3  and  4.92  p.m2  pun-3,  respectively.  These 
values  are  comparable  to  those  reported  by  Wilson  et  al.  [8,9]  who 
calculated  the  total  surface  areas  per  unit  volume  for  a  Ni:YSZ  com¬ 
posite  electrode  as  3.0  pun2  p,m-3  for  YSZ,  1 .6  p,m2  p.m-3  for  Ni  and 
2.4  p.m2  p.m-3  for  the  pore  space. 

3.3.  Effective  transport  properties 

Two  computational  domains  of  two  different  sizes  are 
selected  from  the  3D  reconstructed  volume  of  the  LSM  and 
pore  spaces.  The  3D  computational  domains  have  the  sizes 
of  9.12  pun  x  4.75  pun  x  4.08  p.m  denoted  as  ‘Sample  A’  and  of 
4.98  p,m  x  3.04  p,m  x  4.08  pun  denoted  as  ‘Sample  B\  Fig.  13  shows 
the  computational  domains  required  to  calculate  the  transport 
properties.  Fig.  13(a)  and  (b)  displays  the  computational  domains 
for  Sample  A  and  Sample  B,  respectively.  The  3D  computational 
domains  shown  in  Fig.  13  were  constructed  by  the  finite  vol¬ 
ume  cut-cell  approach  through  MicroFOAM.  To  evaluate  effective 
transport  properties,  we  consider  the  present  FIB-SEM  recon- 
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Table  1 

Normalized  effective  properties  in  x,  y,  and  z  directions  based  on  the  FIB-SEM  recon¬ 
struction  and  the  numerical  reconstruction. 


Sample  A 

Sample  B 

FIB-SEM 

Numerical 

FIB-SEM 

Numerical 

kf/k 

0.242 

0.159 

0.265 

0.184 

kf/k 

0.225 

0.176 

0.294 

0.207 

kf/k 

0.146 

0.174 

0.165 

0.214 

Df/D 

0.304 

0.298 

0.276 

0.265 

Df/D 

0.297 

0.305 

0.284 

0.270 

Df/D 

0.225 

0.312 

0.211 

0.290 

structed  images  as  well  as  the  numerically  reconstructed  structures 
reported  in  Kenney  et  al.  [3].  The  numerically  reconstructed  elec¬ 
trode  assumes  randomly  distributed  LSM  particles  laid  out  to 
achieve  a  target  porosity  of  50%.  Both  reconstructed  volumes  are 
converted  into  the  body-fitted  cut-cell  based  finite  volume  meshes 
using  the  MicroFOAM  [12].  The  normalized  effective  electronic  con¬ 
ductivity  and  the  normalized  effective  gas  diffusivity  are  calculated 
by  Eq.  (5)  using  both  reconstruction  techniques.  Table  1  summa¬ 
rizes  the  normalized  effective  properties  in  x ,  y ,  and  z  directions 
calculated  from  two  different  sizes  of  computational  domains. 

Several  noteworthy  results  can  be  observed  from  Table  1.  Both 
the  electronic  conductivity  and  the  gas  diffusivity  are  calculated  for 
each  direction  of  the  phase  domains.  It  is  observed  that  the  effec¬ 
tive  properties  calculated  by  FIB-SEM  reconstruction  seem  to  be 
direction  dependent  contrary  to  those  evaluated  by  the  numerical 
reconstruction.  For  the  FIB-SEM  reconstruction  case,  the  normal¬ 
ized  effective  properties  in  z-direction,  i.e.,  kff/k  and  Dfff/D,  are 
roughly  25-44%  smaller  than  the  corresponding  properties  in  x  and 
y  directions.  This  difference  is  significant  and  reveals  the  anisotropy 
of  the  measured  sample.  In  comparison,  the  results  for  the  numeri¬ 
cally  reconstructed  samples  vary  by  less  than  1 0%.  The  anisotropy  of 
the  measured  sample  may  have  some  bearing  on  the  milling  direc¬ 
tion  at  the  time  of  performing  FIB-SEM.  It  is  also  likely  due  to  the 
presence  of  a  large  crater,  as  observed  in  Fig.  12.  It  is  to  be  noted 
that  such  a  local  feature  is  not  taken  into  account  in  the  numer¬ 
ically  reconstructed  samples.  Flence,  these  observations  question 
the  ability  of  the  numerical  models  to  represent  the  actual  structure 
of  an  electrode  adequately.  Flowever,  one  should  keep  in  mind  that 
the  numerically  reconstructed  samples  are  an  average  description 
of  the  electrode  microstructure.  Kenney  et  al.  [3]  recommended  the 
use  of  at  least  20  numerical  structures  in  order  to  obtain  converged 
values  of  the  microstructural  properties.  Similarly,  more  samples 
would  be  necessary  over  which  a  series  of  FIB-SEM  images  are  to 
be  conducted  in  order  to  obtain  properties  that  do  not  depend  on 
local  features  such  as  cracks  or  large  craters.  However,  the  latter 
proposition  is  time  consuming  and  often  prohibitively  expensive. 

For  both  the  actual  sample  reconstructed  by  FIB-SEM  and  the 
one  created  numerically,  the  effective  properties  exhibit  larger  vari¬ 
ations  in  the  three  directions  when  the  sample  size  is  small.  The 
importance  of  having  a  sufficiently  large  sample  size  was  already 
acknowledged  in  previous  experimental  works  by  Wilson  et  al.  [30], 
Iwaietal.  [ll],Gostovicet  al.  [31].  Wilson  et  al.  [30]  recommends  to 
use  a  volume  of  at  least  500  |xm3  for  the  sample  to  yield  microstruc¬ 
tural  properties  that  are  representative  of  the  overall  electrode. 
Additionally,  Kenney  et  al.  [3]  and  Metcalfe  et  al.  [14]  found  that 
the  size  of  the  domain  should  be  at  least  14d  or  15d  in  each  direc¬ 
tion,  where  d  is  the  mean  particle  diameter,  for  the  calculation 
of  the  triple-phase  boundary  length.  For  the  electrode  considered 
here,  this  would  mean  a  domain  volume  greater  than  1000  [xm3. 
Choi  et  al.  [12]  found  that  a  size  of  lOd  x  lOd  x  lOd,  or  390  |xm3  in 
our  case,  is  sufficient  for  the  calculation  of  the  effective  properties. 
Therefore,  it  can  be  commented  that  the  volume  of  the  sample  mea¬ 
sured  by  FIB-SEM  (approximately  177  |xm3)  in  the  present  study 


may  be  smaller  than  required  based  on  the  above  arguments.  One 
needs  to  be  cognizant  of  such  volume  dependencies  on  transport 
properties  calculated  based  on  FIB-SEM  images. 

4.  Conclusions 

In  this  work,  a  dual-beam  FIB-SEM  technique  was  applied  to 
reconstruct  the  microstructure  from  a  porous,  ceramic  LSM  cathode 
electrode  material  for  SOFC.  Preparation  of  LSM  sample  and  acqui¬ 
sition  of  2D  cross-sectional  images  from  the  FIB-SEM  facility  was 
provided.  Then,  a  formulated  methodological  approach  was  used  to 
process  the  FIB-SEM  2D  images  for  reconstructing  3D  microstruc¬ 
tures.  For  the  first  time,  it  was  found  that  the  image  post-processing 
could  substantially  affect  the  geometric  properties  such  as  porosity 
s  and  internal  surface  area  S.  The  choice  of  threshold  values  as  well 
as  median  filter  operations  affected  variations  of  the  porosity  and 
internal  surface  area  up  to  33%  and  25%,  respectively.  It  is  observed 
that  not  applying  the  MO  operation  significantly  overestimates 
the  porosity  (by  about  10%  for  r  =  2  pixel)  and  the  internal  surface 
area  (by  about  25%  for  r  =  2  pixel).  The  resolution  of  the  surface 
triangulation  used  in  3D  reconstructions  also  varied  the  poros¬ 
ity  s  and  internal  surface  area  5  by  14.5%  and  4.4%,  respectively. 
These  variations  indicate  the  sensitivity  of  imaging  processing  and 
phase  segmentation  methods  on  calculating  different  properties. 
3D  tessellation  is  also  important  and  must  be  done  with  adequate 
resolution. 

The  resulting  3D  reconstructed  volumes  were  converted  into 
the  computational  domains  based  on  the  body-fitted  cut-cell  based 
finite  volume  meshes.  The  normalized  effective  transport  prop¬ 
erties  in  x,  y,  and  z  directions  were  calculated  on  computational 
domains  reconstructed  by  the  FIB-SEM  technique  as  well  as  by  the 
numerical  model.  It  was  observed  that  the  effective  properties  cal¬ 
culated  by  FIB-SEM  reconstruction  are  more  anisotropic  than  those 
evaluated  by  the  numerical  reconstruction.  This  local  anisotropy  of 
the  measured  sample  was  likely  due  to  the  presence  of  the  large 
crater  (as  shown  in  Fig.  12)  as  well  as  a  smaller  sample  volume 
used  in  the  current  FIB-SEM  reconstruction.  It  is  suggested  that  a 
large  sample  volume  and  a  statistically  large  number  of  samples 
may  remove  such  local  anisotropy  behaviour. 
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